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Abstract 

The production of (6c)-quarlconium {B(. meson and its excited states) or (c6)-quarkonium [B^ 
meson and its excited states) via top quark t or top anti-quark i decays, t — > (6c) + c + or 
i (c6) + c + W~ , respectively is studied within the framework of NRQCD. In addition to the 
production of the two color-singlet S-wave states |(6c)(^S'o)i) or | (c6)(^S'o)i) and |(6c)(^5i)i) or 
\{cb){'^Si)i), the production of the P-wave excited (be) or (c6) states, i.e. the four color-singlet 
P-wave states \{bc){^Pi)i) or \{cb){^Pi)i), and {bc){^Pj)i) or {cb){^Pj)i) (with J = (1,2,3)), is 
also studied. According to the velocity scaling rule of NRQCD, for the production of P-wave 
excited states the contributions from the two color-octet components \{bc)(^ So)s) or \{cb)(^ So)s) 
and \{bc){^Si)s) or |(c6)('^S'i)8) are also taken into account. We quantitatively discuss the 
possibility and the advantages in experimental studies of Be or Be meson and its excited states via 
the indirect production at LHC in high luminosity runs and at LHC possible upgraded versions 
such as SLHC, DLHC, TLHC etc. in future. 
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I. INTRODUCTION 

Since Be meson was observed by CDF in 1998 [ll, the interests in studying of the meson 
Be are increasing and a new stage of the studies has been started. As pointed by the authors 
of Refs. , the product cross-section at Large Hadron Colhder (LHC), CERN, is much 
greater than that at Tevatron, Fermilab, therefore, more precise experimental studies of 
Be meson are expected at the forthcoming of LHC. Indeed, some progress at LHC in the 
experimental study of Be physics can be achieved, especially at the beginning stage of LHC 
(LHC just runs at 'low' luminosity L = 10^^ ~ 10^^cm~^s~^). However, when LHC runs 
at higher luminosity, such as the design luminosity L = lO^^cm'^s^^ later on, due to the 
requests coming from the main purposes of LHC (such as searching for Higgs particle and 
SUSY partners etc) and the restricts from the abilities to record the events for the detector 
(e.g. the detector cannot record too frequently events etc), the condition for triggering the 
events occurring in collisions has to be set such that too many Be events via the direct 
hadronic production according to the theoretical estimate of 'direct' i?c-productionj3, sl, I4, 
are cut off. As a result, the direct hadronic production of Be cannot be expected to 
make much progress in i?c-meson study in the high luminosity stage of LHC runs. Baring 
the situation pointed out here and the possible upgrade for LHC (SLHC, DLHC and TLHC 
etcj3]) in mind, the possibility to study Be meson experimentally via indirect production 
of Be meson, namely, via producing a huge amount of top antiquarks t and their decays, 
is worthwhile to think seriously about. It is because that at LHC no matter how high the 
luminosity, the produced top quark(s) shall never be cut off by the trigger condition set 
down for any experimental purposes, and the frequency of the t-quark production can be 
stood up for the detector always. Furthermore, the mechanisms for producing the doubly 
heavy mesons, such as rje, rji,, Be, ■ ■ ■ and J/ip, T, B*, ■ ■ -, are interesting too. The indirect 
production of Be or Be{oT B~), including its excited states, via t-decays or t-decays may 
offer some knowledge on the mechanisms, therefore, this paper is devoted to study the 
indirect production of Be or Be meson via t-decays or t-decays. Without confusing and 
for simplifying the statements, later on we will not distinguish Be and Be unless necessary, 
and all results for Be and Be obtained in the paper are symmetric in the interchange from 
particle to anti-particle. 

The doubly heavy meson production via top quark decays is special interesting from the 
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point view of precise testing perturbative quantum chromodynamics (pQCD) [8|] . The meson 
Be is the ground state of the heavy-flavored binding system (6c), and it is unique 'doubly 
heavy-flavored' anti-meson in Standard Model and is stable for strong and electromagnetic 



interactions. Although in literature the 'direct 



ladornic production of Be meson has been 
(references therein) and CDF discovered 



thoroughly studied, e.g. see Refs. 0, y, Q, y, B 

the meson which just come from the 'direct' production, as a compensation to understand 
the production mechanisms, it is quite interesting that to study the production of Be meson 
'indirectly' through t-quark decays, especially, considering the fact that numerous t-quarks 
may be produced at LHC. The theoretical studies of the direct production [2], y, Ij, 5, 6] is 
based on NRQCD j^, so now we study the indirect production based on NRQCD too. 

In the framework of effective theory of NRQCD, a doubly heavy meson is considered as 
an expansion of a series Fock states. The relative importance among the infinite ingredients 
is accounted by the velocity scaling rule. Namely the physical state of B^ B*, hj^^ and 
can be decomposed into a series of Fock states as follows: 

\Be) = 0{v')\{bc),CSo)) + 0{v')\{bc)sCPi)9) + ■■■ 

m = 0{v')\{b-cU'S,)) + 0{v')\{b-c)sePj)g) + ■■■ (1) 

and 

\hsj = Oiv'')\{bc)^CP^)) + Oiv')\ibc)sCSo)g) + ■■■ 
\xi) = 0{v')\{b-chePj)) + 0{v')\{b-c)seS,)g) + ■ • • , (2) 

here v is the relative velocity between the components. The thickened subscripts of the 
(be) denote for color indices, 1 for color singlet and 8 for color-octet; the relevant angular 
momentum quantum numbers are shown in the parentheses accordingly. According to the 
velocity scaling rule of NRQCD, the probability of each Fock state in the expansion is 
proportional to a definite power in v as indicated as that in Eqs. f|T][2|) . Since the value 
of is around 0.1 ~ 0.3 [9|], which is not too small, and the contributions from the two 
color-octet S-wave components to the P-wave production might be comparable with those 
from the color-singlet components. So we shall consider the two color-octet components: 
\{bc)sC'So)g) and \{bc)s(^Si)g) , in addition to those color-singlet components in the mesons 
Be, B*, and xi- 
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The calculations of the process are very complicated and lengthy by using the conventional 
trace techniques to calculate the amplitude square due to the two massive particles and 
bound state effects. To shorten the calculations and to make the results more compact, we 



adopt the method used in Ref. lOj] to do the calculations. For convenience, we will call it 
as the 'new trace amplitude approach'. Under this approach, we first arrange the whole 
amplitude into several orthogonal sub-amplitudes Mss' according to the spins of the t-quark 
(s') and c-quark (s), and then do the trace of the Dirac 7 matrix strings at the amplitude 
level, which result in explicit series over some independent Lorentz-structures, and finally, 
we obtain the square of the amplitude. During the calculating, some useful tricks have also 
been introduced to make the expressions more compact. More detail of the techniques and 
all the necessary expressions for the amplitudes of t(po) ipc){pi) + 0(^2) + W^ij)^) with 
(fec)-quarkoniumin |(6c)(i5o)i), \{hc){^ Si)^) , |(6c)(iPi)i) and {hc){^Pj)i) (with J = (1,2,3)) 
respectively are put in the APPENDIX B. While the expressions for the contributions from 
the two color-octet components | (6c) (^6*0)8) and | (6c) (^5*1)8) can be obtained from those 
from the two color singlet S-wave components by changing the overall color-factor and the 
corresponding matrix elements. This approach is different from that of the spinor techniques 



or the so called 'helicity amplitude approach', which has been proposed in Ref.[llJ and 
improved in Ref.|l2|]. Under the 'helicity amplitude approach', one can also derive compact 
results that can be further expressed by spinor-products at the amplitude level and then 
do the numerical calculations ^. The above two methods are complement to each other 
and both can derive simple and compact expressions at the amplitude level. The 'helicity 
amplitude approach' is more suitable for the numerical calculations and is more quicker 
for calculations, since full components of the helicity amplitude can be evaluated at the 
amplitude level. While for the 'new trace amplitude approach', at the amplitude level only 
the coefficients of the basic Lorentz structures are numerical. However from the amplitude 
derived from the 'new trace amplitude approach', one can sequentially result in the squared 
amplitude, which is more easier to be compared with the results derived by the traditional 
trace techniques. 



According to Refs.|13l|. one may expect at LHC to produce ~ 10 tt-pairs per year 



^ Here, we refer to Ref. 5] for an example on the 'helicity amplitude approach', where full processes of the 
approach from the formulae deduction to the numerical calculation can be seen explicitly. More over some 
tricks to simplify the massive amplitudes can be found there. 
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under the luminosity L = lO^^cm ^. Considering the possible upgrade for LHC and the 



estimate by Refs. 



131 ]. we will assume that one may obtain ~ 10^ — 10^° tt-pairs per year to 



examine the possibility to observe the meson Be via decay of the produced t t-pairs precisely 
and to examine the advantages in the indirect way to observe the meson and its excited 
states. 

The paper is organized as follows. In Sec. II, we show our calculation techniques for the 
process t(po) {bc){pi)+c{p2)+W^{p3). Then we present numerical results and make some 
discussions on the properties of the (6c) -pro duct ion through t-quark decays in Sec. III. The 
fourth section is reserved for a summary. All necessary expressions are put in the appendices 
finally. 



II. CALCULATION TECHNIQUES 

Under the NRQCD framework [o], [l^, the total decay width for the production of {be 
quarkonium through the channel t{po) (6c) (pi) + c(p2) + W~^{p3) takes the form: 

, {On) 



T = J2HrXt^ibc) + c + W' 



N, 



col 



(3) 



where Ncoi refers to the number of colors, n stands for the involved state of 6c-quarkonium. 
Ncoi = 1 for singlets and Ncoi = — 1 for octets. With the help of the saturation approxi- 



mation j^, the involved decay matrix elements can be written as ^, 

1 



{bcCSo)i\OiCSo)\bcCSo)i) 
{bcCSi)i\OiCS,)\bcCSi)i) 
{bcCP,),\OiCP,)\bcCP,),) 
{bci'Po)i\0,CPo)\bcCPo)i) 
(6c(3Pi)i|0(3Pi)i|6c(=^Pi)i) 

(6c(3p2)i|Oi('P2)|6c(3P2)i) 



1 



\xt^,\bcCSo)mi + 0{v% 
\xt^A\bceSi)imi + 0{v% 



(4) 
(5) 
(6) 



^ {0\xti-'7;^-^)A\bci'Po)i)ni + 0{v% (7) 



1 



2 

i . 



\Xc 



Bxa)i;,\bc^CP^),)\'[l + Oiv% (8) 



^ X^(-iSWa(^-))V^,|6c(3p2)i)p[l + Oiv% (9) 



where the subscript 1 or 8 indicates that the operator is a color singlet or a color octet, 
ip~^ is the Pauli-spinor field that create a heavy quark, x is the Pauli-spinor that create 



Here as suggested in Ref. 1J|, an overall factor l/2Nc is introduced into the color-singlet matrix elements. 
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a heavy antiquark, — + igA^ is the gauge-covariant derivative, A is the SU{3)- 
matrix- valued gauge field. The operator D is the difference between the covariant derivative 
acting on the spinor to the right and on the spinor to the left, which is defined by x '^'^ = 
X^(Dip) — (DxV^- The notation T^^^^ is for the symmetric traceless component of a tensor: 
rpiij) ^ (jij ^ yii)/2 - T''^5'K Furthermore, we have 

^ -\xti^,\h-cQS,),) = -4^Rsim + 0{v')] (10) 



^ mxtcJ^I^^M'SM^)) = ^i?5(A)e[l + 0{v^)\ (11) 



^ -MxM-i-m^m'PMe)) - \/5^p(A)4i + o{v')] (12) 



y2iV,^ i/vc V 2 ^ ^^^^ ^' V 



^ ^ :(0|Xc^<t(-^D X <7)^,|6c(^Pi)i(6)) = ./5i?^(A)6[l + 0(^^)] (14) 



V2 V2iV: '-^'^ ' 2 ' '''' " V47r 

^(0|Xc^<7(^B«<7(^))V'.|6c(3P2)i(6)) = y£i?'p(A)6^^[l + 0(^2)], (15) 

where -R5(A) is the average radial wavefunction for \S state averaged over a region of size 
1/A centered at origin, ^p(A) is the average derivative of the radial wavefunction at origin 
of size 1/A. For convenience, we shall take ^5 (A) and ^p(A) to be the phenomenological 
values -R5(0) and -Rp(O), which may be derived from the QCD potential models and relate 
to certain observable such as the width for electromagnetic annihilation etc. 

Although we do not know the exact values of the two decay color-octet matrix elements, 
(6c("'^5'o)8|C^8(^'S'o)|&c(^S'o)8) and (^ci^ S\)^\0%i^ S\)'^c{^ S\)%) ^ we know that they are one order 
in higher than the 5'-wave color-singlet matrix elements according to NRQCD scale rule. 
More specifically, based on the velocity scale rule, we have: 

(6c(i^o)8|08C^o)|fc(i^o)8) ^ A5(t;)^-(6c(^5o)i|0i(%)|?>c(%)i) 
{h-ci^SMO^i^S^)\hci^S^\) ~ A5(^;)' • {h-ci^SMOi^S^)\hci^S^y) , (16) 

where the second equation comes from the vacuum-saturation approximation. A5(t') is of 
the order or so, and we take it to be within the region of 0.10 to 0.30, which is in consistent 
with the identification: iS.s{v) ~ agiMv) and has covered the possible variation due to the 
different ways to obtain the wave functions at the origin (5"- wave) and the first derivative of 
the wave functions at the origin (P-wave) etc. 
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FIG. 1: Feynman diagrams for t{po) {bc){pi) + c{p2) + W'^{p3), where 6c-quarkonium is in eight 
states, i.e. the six color-singlet states (^^o)!, ('^5'i)i, {^Pi)i and (^-Pj)! (with J = (1,2,3)), and 
two color-octet states (""^ 6*0)8 and {'^Si)s respectively.. The left and the right diagram are for the 
hard scattering amplitudes ^1 and A2 respectively. 

The Feynman diagrams for the process of t{po) — >■ bc{pi) + c{p2) + W~^{p3) are shown in 



Fig.([T]), where 6c-quarkonium is in eight states: six color-singlet components {^S{ 



01, 



(}Pi)i and (^Pj)i (with J = (1,2,3)), and two color-octet components (^6*0)8 and ('^S'i)8 
respectively. Based on the phase-space integration simplification as shown in the Appendix 
A, the decay width of the process can be written in the form: 

3|Mp 



dT 



2567r^m? 



X — asias2, 



N, 



col 



(17) 



where the extra factor 3 in the numerator comes from the sum of the c-quark color, |Mp 
is the mean square of the hard scattering amplitude, i.e. |Mp = S l-^i + -^2^ with Ai 
and A2 are two amplitudes of the process, si = {pi + P2Y and S2 = {p2 + PzY- 

We deal with either L = or L = 1 state here. The two hard scattering amplitudes of 
the process, which correspond to the left and right diagram of Fig.([l]), can be written as 



A 



5=0,L=0 



A^ 



5=0,L=0 



and 



S=1,L=0 



A 



S=l,L=a 



-- iCui{p2,s) 
-- iCui{p2,s) 

iCui{p2,s) 
iCui{p2, s) 



7m 



KM 



{p2 + puy'^^ {pi + - ml 



Apz)Pl 



(P2 + Pur {pi2 + ps) - l^t 

£^(pi)n-(g)_ + 



2 7m 



{P2+Pliy '^^{Pl+P2 

g^(pi)n° (g) 
7m nP3)PL 



mt 



Ap^)Pl 



'12+ p-i + mt 



{P2+Piiy 



{Pi2+PzY - m 



2 7m 



Mj(pO,s')|g=0 (18) 

Uj{po,s')\q=Q (19) 
%(po,s')|g=o (20) 

Mj(pO,s')|g=0 (21) 



and 



5=0,L=1 



iCe1{pi)ui{p2,s] 



d 
dqa 



7m 



h+ P2 + mb 



KM _ 

{P2 + Pu^'^'^'iPl + P2) 



mt 



Ap^)Pl 



Uj{po,s%i^2) 



and 

i5=l,L=l 



dqa 



{P2 +Pu)^ 



AT 



iCeip{px)ui{p2,s) 
iCeip{pi)ui{p2,s) 



d 

d 
dqa 



7m 



'(pi2 +P3)^ - mf 



h+ P2 + rrib 



Uj{po,s%2d) 



[P2+Piiy {Pi+P2)~rnf, 



(q) 

{P2+PllY 



'12+ ^3 + m 



{pu+Psf 



i7m 



ujipo,s')\iM) 
Uj{po,s'){^) 



with the color factor C = or Co for color-singlet and color-octet respectively, Cs — j^^ij 
and Co = ^(\/2T»T^r«)ij (\/2T^ stands for the color of the color-octet be state), ^(pa) is 
the polarization vector of W^, Pl = ^—2^ and Pr = q, pii and pi2 are the relative 

momentum between the two constitute quarks of (6c)-quarkonium and the momenta of these 
two constitute quarks respectively. More explicitly, we have 



Pii = ^Pi + q and P12 = —P2 - q, 



(26) 



where M ~ nih + rric- s^ipi) and e"(pi) are the polarization vectors relating to the spin 
and the orbit angular momentum of (6c)-quarkonium, ^^^(pi) is the polarization tensor for 
the spin triplet P-wave states with J = 0, 1 and 2 respectively. The covariant form of the 
projectors can be conveniently written as 

fM 



KM)- 



{Ai - "^0)75(^12 + "^6), 



and 



KM) = 7:^(^11 - ^c)7"(^i2 + m,), 



To do the simplification of the projector, the following simplification shall be useful: 



n°(o) 



1 



2VM 



75(^1 + M) , n°(o) 



1 



2VM 



7"(^i + M), 



and 



4-nt 



M 



7a 7/3(1^1 + mb- rric) - 2gap(fn - "f^c) 



(27) 



(28) 



(29) 



(30) 



(31) 



dqa^ Pl -i/iy-" Arribmc 
Here the properties: p"£" = and p"£"^ = pfs"'^ = 0. pJi — jiPi P12 — 1a P^^ are 
applied. After substituting all these relations into the amplitudes and doing the possible 
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simplifications, the amphtudes then be squared, summed over the freedoms in final state and 
averaged over the ones in initial state. And the selection of the appropriate total angular 
momentum quantum number is done by performing the proper polarization sum. If defining 

n«/3 = -gap + , (32) 

the sum over polarization for a spin triplet S-state or a spin singlet P-state is given by 

^^f^a^Q,' = nQ,Q,/, (33) 

where = or 1^ respectively. In the case of ^Pj states, as for the three multiplets £^a^(pi) 
with J = 0, 1 and 2, the sum over the polarization is given by 



1, 

-(1)^(1)* 

2' 



^al^a'P' - -'^af3^a'f3' (34) 
1 



Xl^a/S^a'/?*' - ;j(naa'n/3^/ - Ea/j/IlQ,/^) (35) 
1 1 

All the terms for the squared hard scattering amplitudes |Mp for S-wave and P-wave 
states are put in Appendix B accordingly, where the detail of the calculation techniques are 
also attached for convenience. To effectuate the calculations and to make the results more 
compact, we adopt the method used in Ref . [lol| . As stated in the Introduction, we call it 
the 'new trace amplitude approach' for convenience. Under the approach, we arrange the 
whole amplitude into several orthogonal sub-amplitudes Mgg' according to the spins of the 
t-quark (s') and c-quark (s) first, and then do the trace of the Dirac 7 matrix strings at the 
amplitude level by properly dealing with the massive spinors, which result in explicit series 
over some independent Lorentz-structures. The expressions for these coefficients of all the 
considered channels are put in Appendix B. With the help, one can do the square of the 
amplitude easily. As a cross-check of our results, we adopt the traditional trace techniques 
and also the FDc[l^ package to derive the numerical results of the mentioned processes. 
Indeed we find a well agreement among these methods. 

As a comparison and for later usages, let us present the width for the two body decay 
t{pi) b{jp2) + W^ij)^), which is dominant for the t-quark decay: 

\^i_y2Y + :,^^i + y-^2x') , (37) 



4V27r 

where IP2I = 0- ~ — — {x + y)^), rriu, = rritx and rrii, = rrity. 



III. NUMERICAL RESULTS OF DIFFERENTIAL CROSS-SECTIONS 



In numerical calculations, we take the parameters as follows: 

nib = 4.9 GeV, rric = 1.5 GeV, mt = 176 GeV, = 80.22 GeV, a,(2mc) = 0.26, (38) 



{be 


) + W+ + 


c are: 




.(6c)[(35i)i] 


= 0.79 MeV 






= 0.57 MeV 


r*^ 


>(6c)[(iPi)i] 


= 0.057 MeV 


r*^ 




= 0.034 MeV 


r*^ 


.(6c)[(3Pi)i] 


= 0.070 MeV 


r*^ 


>(6c)[(3P2)l] 


= 0.075 MeV 


r*- 




= 0.091 X -y^ MeV 




■(6c)[(iSo)8] 


= 0.070 X MeV 



and g — 2^/2mw^J Gp/ \/^- Then, the decay width oit^ W'^ + & is 

T{t ^ + 6) = 1.59GeV. (39) 

(40) 
(41) 
(42) 
(43) 
(44) 
(45) 
(46) 
(47) 

where ~ (0.1 ~ 0.3). It may be more useful to show the ratio between the decay width 
of t — (fee) + + c and t — > + 6, since uncertanties from the electro-weak coupling 
can be cancelled out: 

4.97 X 10"^ (48) 
3.58 X 10-^ (49) 
0.36 X 10"^ (50) 
0.21 X 10"^ (51) 
0.44 X 10-^ (52) 
0.47 X 10-^ (53) 
0.57 X 10-V (54) 
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FIG. 2: The invariant mass distributions for t(po) ~^ ibc){pi) + c{p2) + W~^{p3). The left is for 
dT/dsi and the right is for dT/ds2- The diamond hne, the cross hne, the dotted hne, the sohd 
hne, the dashed hne and the dash-dot hne are for (6c)-quarkonium in Fock states: (^Si)!, (^S'o)i, 
(^^2)1, (^A)i, ('i^i)i and (3Po)i respectively. 
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FIG. 3: The differential distributions dT/dcos6i3 (Left) and dT/dcos6i2 (Right) for t(po) — ^ 
(6c) (pi) + c{p2) + W^{p3). The solid line, the dotted line, the diamond line, the dashed line, the 
triangle line and the dash-dot line are for (6c)-quarkonium in Fock states: ('^5'i)i, (^5o)i, (^P2)i, 
(^Pi)i, {^Pi)i and (^Po)i respectively. 



t^(bcWSo)s] 



0.44 X IQ-V. 



(55) 



T{t -^W+ + b) 

Let us show some more characteristics of the decay t —>■ {be) + + c. The differential 
distributions of the invariant masses Si and S2, i.e. dT/dsi and dT/ds2 are shown in Figj2l 
While the differential distributions of cos^ia and cos ^12, i.e. c/F/fi cos 6'i3 and dV /d cos 9\2 is 
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shown in FigjSl where ^13 is the angle between pi and ps, and 612 is the angle between pi and 
P2 respectively in the t-quark rest frame (po = 0). It can be found that the largest differential 
cross-section of dT /d cos 613 is achieved when 613 = 180°, i.e. the (6c)-quarkonium and 
moving back to back in the rest frame of t-quark. And the largest differential cross-section 
of dr/dcos6i2 is achieved when 612 = 0°, i.e. the (6c)-quarkonium and c-quark moving in 
the same direction. 

IV. DISCUSSIONS AND SUMMARY 

In the present paper, we have studied the decay channel t(po) ~^ bc{pi) + 0(^2) + W~^{p3) 
in the leading as calculation but with the f ^-expansion up to v'^, where 6c-quarkonium is 
in one of the eight states: the six color-singlet states (""^ 5*0)1, ('^ 5*1)1, (}Pi)i and {^Pj)i (with 
J = (1,2,3)), and two color-octet states (^5*0)3 and (^5*1)3 respectively. In literature, only 
(^5*1)1 state has been studied Q], however it can be found that all the other considered states 
can be sizable in addition to the (^5i)i state. 

As mentioned in the Introduction, about 10^ t i per year will be produced in the stage of 
high luminosity run at LHC and t-quark events always trigger the detector, then according 
to the present estimate it is possible to accumulate about 10^ events a year via t-quark 
decay at LHC. Moreover, the indirect production of the (6c)-quarkonium may be traced back 
to t-quark decay and has the characteristics in ^^13 , 612 shown in Figs. (2,3) etc, that may be 
used to identify the (6c)-quarkonium events. Thus there may be some advantages in (6c)- 
quarkonium studies via the indirect production in comparison with the direct production. 
Especially when LHC is really upgraded to SLHC, DLHC and TLHC, so 10^ ~ 10^° tt per 
year may be produced, one may expect much progress in (6c)-quarkonium studies is achieved 
then. 

Acknowledgments: This work was supported in part by the Natural Science 
Foundation of China (NSFC) and by the Grant from Chongqing University. 
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APPENDIX A: FORMULAE FOR THE PHASE SPACE INTEGRATION 



In this section, we derive the phase space of i(mt,po) {bc){mi,pi) + c(m2,P2) + 
(7713,^3). The decay width is proportional to the phase space: 

^rocignp||^MV(p„-p.-p,-P3). (Ai) 

where Pi — {Pi,Pi) — (p'i , p] , P^ , P^) ■ Furthermore, in the rest frame of top quark (pq = rrit), 
we have 

« .0 In ,, dWp2d'p3S{pl - ml)S{pl - ml)S{pl - ml)9{p'M) ■ 



dsids2 (2mt)(27r)5 

0{pI)5{Si - {pi +P2)^)5(S2 - iP2+P3y)S'^iP0-Pl -P2 ' Ps) 

oc (2mtX27r)5 ^^^^^^^^'^^^^ ~ ^i)^((Po - Pi - P^f - ml)5{pl - ml)e{p\) ■ 
9{pl)9{mt -pI- pI)5{s, - (po - Pzf)Ks2 - (po - Pi)') 

« ^^^.d'prd'mpf - p1 - ml)8{pf - p1 - ml)e{p\)e{pl) ■ 
9{mt -Pi- pI)S{s2 + ml-ml- 2mtpl + 2pSp3 - 2pi • Ps) 

(X ^^^,dn, sin 9,3d9,sdh9{p'i)9{pl)9{mt - p? - p°) • 

5{si + S2- ml -ml + 2p?p° - 2|pi| • Ips] cos ^13) (A2) 
^^^(^^?)^(^^3)^K -P? -P3)^(^), (A3) 



where p\ — and pg = "^*~2^ ■> Ip^I ~ \Ip^ ~ |Pil — \l P^^ ~ ^i- The step 

function ^(X) is determined by ensuring | cos ^13] < 1, where 

si + S2 - mf - ml + 2p?po 
cos ^13 = oi-ii-i ■ 

2|Pl||P3| 

All these step functions leads to the integration ranges: 



S2-ml + m?)(s2 - "^i + m|) + ^A(s2, m?, m?) A(s2, mi, mi) 



= + _ 1^ ^ — — ^A4) 



= ml + ml 



2s 

[s2 -ml + ml){s2 -ml + m|) - ^A(s2, ml, mf)A(s2, ml, ml) ^ 



sf^ = (m2 + m3)' (A6) 
s^^^ = [mt-mif (A7) 

where A(x, y,z) = {x — y — zf" — Ayz. 
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Furthermore, we can obtain the cos ^13 distribution from Eg. (1A2P : 
where the extra Jacoiban 

-\pi\\n\ 



J 



1 _ P-j J \p3.\{rn{^ +mt'^ ~ S2) cos 013 



^'"i mt -\/mi4+{mt2-s2)^-2mi2(jnt2+s2) 



(A9) 



and 



IpsI cos 6*^3 - {mt - pi) 



mtpl) - mi^pl) - mi|p3| cos6'i3 - 2mi\ni2^ + 2{mt - plf - Si - 2|p3|^ cos6'i3) 



1/2 



(AlO) 



The 9{X) function determines the boundary of si. 

SVmin <si< {nit - m^if (All) 

where 



with 



m2^mt^ + mi'{mt^ cos 6\.^ + m3^(cos 6'^3 — 1)) + mimtvY t 
mf — mi (^i — cosfc'i3j 



y = {rrii^ — VTi'i' + m3^ — m^f' — [vtix^ — 1{m\ — + mf)m1 + 

- ^3)^ - mMm. + m,f - ml)) cosOl, + Am.^m,^ cosOt,. (A13) 



It should be noted that the integration over cos0i3 should be from 1 to —1. The distribution 
for $12 can be obtained in a similar way. 

APPENDIX B: AMPLITUDE OF THE PROCESS t{po) {bc){pi) + c{p2) + VF+(p3) 

Now let us illustrate the method to calculate the amplitude squared. In general, the 
amplitude for the process in which two massive fermions with spin projection s and s' are 
involved can be written as 

Mss' = Us{P2)AUs'{Pq), (B1) 
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where p2 and po denote the momenta of the top quark and the charm quark with mass rrit 
and and A is an exphcit string of Dirac 7 matrices for the process, which can be read 

Let us introduce a massless spinor m^(/co) with a hght-hke momentum and negative 
hehcity first. Thus U-{ko) is satisfied with the following projection relation: 

m_(A;o)m-(A;o) = tu- ^0, (B2) 

where U- = (1 — 75) /2. By introducing another spacelike vector ki which satisfies the 
relations: 

fci ■ fci = -l,fco ■ fci = 0, (B3) 
then the other, massless, independent and positive helicity spinor u+{ko) may be constructed: 

u+{ko) =^iM_(fco). (B4) 
It is easy to check that u^{ko) is satisfied with the projection 

u+{ko)u+{ko) = UJ+ ^0, (B5) 

where a;+ = (l + 75)/2. Using these massless spinors, one can construct the massive spinors 
for the fermion and antifermion as follows: 



Usiq) = {^+m)u^{ko)/\/2ko ■ q, (B6) 
u~s{q) = {^+m)u+{ko)/^j2ko-q, (B7) 

with the spin vector s^: 

_ q^ m 
m q ■ kQ 

Using the above identities, we can write down the amplitude M±s±s' with four possible spin 
projections in the trace form for the 7-matrices: 

Mss' = NTr[{^o + mt) ■ 00^ f^o ■ (^2 + m,) ■ A] (B8) 

M_s-s' = NTr[{^o + mt) • ^0 • (^2 + m,) ■ A] (B9) 

M_ss' = iVrr[(A + mt) ■ co^ ^0 h ■ (^2 + m,) ■ A] (BIO) 

Ms-s' = iVTr[(/o + mt) ■u:+hh-{l^2 + m,) ■ A\ (Bll) 
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with the normahzation constant N — 1/ ^ 4(A;o • Po) (^o • P2) ■ Thus the squared unpolarized 
matrix elements can be written as 

|M|2 = \Mss'\^ + \M^s-s'? + + \Ms-s'?. 

Next, we are to simphfy the calculation. For such purpose, we recombine these M±s±s' 



n — 1, 


• • • , 4) as 


follows: 








Ml = 






N 

-j=Tr[{^Q + mt)- 


h ■ (^2 + rUc) ■ A] 


(B12) 


M2 = 




- M_s-S') = 


N 

-j=Tr[{]^Q + mt)- 


hl^ • (^2 + ^c) • A] 


(B13) 


M3 = 






N 


h ¥q ■ (^2 + rric) ■ A\ 


(B14) 


M4 = 






N 


■ 75 h h ■ {i>i + "T.c) • ^] 


(B15) 



and then we have |M|^ — \Mi^ + IM2P + + |M4p. In order to write down A as 

explicitly and simply as possible, we set the vector k^: 

ko=P2- OiPo, (B16) 
where the coefficient a is determined by the requirement that ko be a light-like vector: 

^ ^ P0-P2 _^ jA_ 



with A = \J{pq •P2Y — mf-ml . Furthermore, if choosing ki: ki-po — and ki-p2 — 0. And 
then the resultant M„ can be simplified as 

Ml = LiTr[(^ + mt) ■ (^2 + mc) • A] (B17) 

M2 = -L2Tr[(^ + ■ 75 ■ (^2 + ■ A] (B18) 

M3 = L2rr[(^o + m)- h ■ {i>2 + me) • ^] (B19) 

M4 = Lirr[( A + mt) • 75- ^1 • (^2 + va^ ■ A\ (B20) 



with 



1 1 

2Li = and 2L2 



VPo ■ P2 + mcmt VPo ■ P2 - mcmt 

The value of /ci is arbitrary, and we take its explicit form as 



A;^ = i«e'^'^''>,pipP2a, (B21) 
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where k is a suitable normalization constant and can be expressed as 

h = 1^15 [{po ■ Pi) i)2 + (p2 ■ p\) A - {pi ■ po) i>\- i>T i)v A 

Substituting fci into Eqs. (]B19|B20p . we obtain 



M3 = + fi:[(po ■ Px)m, - (pi ■ P2)mt]M2 (B22) 

M4 = - kHpo ■ pi)m, + (pi ■ p2)mt]Mi (B23) 

where 

= 7^Tr[(A + mt) ■ 75- ^1 ■ (^2 + m^) ■ A] (B24) 

= 7^Tr[(^ + mt)- ■ (^2 + m,) ■ A] (B25) 

The amplitudes for the process t{po) {bc){pi) + 0(^2) + ^^'^(Pa) can be expanded over 
some basic Lorentz structures: 

m m 

M,{n) = Y,A'.{n)B,{n){t = l-4) , M^{n) = J2 {n)B,{n) (z = 3,4) (B26) 

where m is the number of basic Lorentz structure Bj{n), whose value dependents on the 
(6c)-quarkonium state n: e.g. m = 3 for n = {bc)[^So]i, m = 11 for n = (6c)pS'i]i and 
{bc)[^Pi]i and m = 30 for n = {bc)[^Pj]i. As for Aj{n) and A'^{n), they can be expressed by 

A^.{n) = Af{n) + k[{po ■ pi)m, - {pi ■ p2)mt]A^^{n) (B27) 
A^n) = Af{n) - 4{po ■ pi)m, + (pi ■ p2)mt]A]{n) (B28) 

The explicit expression for A^j''^{n) and A^ (n) of each state shall be listed in the following 
subsections. 

To short the notation, we define some dimensionless parameters 

rrib iTLc rriuj M 

n = — , r2 = — , rs = , r4 = — 

rrit rrit rrit mt 



and 



u = pi -p^/m] = ^(^3 -"^1, -ml),v = p2-p3/m'l = ^(^2 - m^-ml), 



2 _ 1 ^_2 , _2 „ \ ^ „ /_2 _ 1 _2 ^2 



2mr - . 2mt^ 

y = Po-Pi/mt = 7^{mt + ml - Sa), 2 = Po ■ ^2/"^* = TT^i^t + - ^3 
2mt 2mt 
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where si = (pi + ^2)^, S2 = (p2 + PsY and S3 = (pi + ps)^, which satisfy the relation: 
si + 52 + S3 = TTif + ml + + mg^. And the short notations for the denominators are 

1 1 ^ _ 1 1 

(P2 + Pll)^ (Pl + P2Y -ml' (p2 + Pll)^ (P3 + Pl2)^ - ' 

mf 1 1 

(p2+P?i)MPi+P2)^-mg ' (p2+P?i)MP3+P?2)^-mr 
1 

(P2+P?l)M(P3+P?2)'-m2)2- 

Furthermore, the following relations are useful to short the expressions: 

y + z + w — 1, X + u + rl — y, x + v + — u + v + r\ — w. 



ds = 
d5 = 



1. Coefficients for spin-singlet S-wave state: (6c)[^S'o]i 

There are three basic Lorentz structures Bj for the case of (&c)[^5'o]i, which are 

^1 = , B2 = , ^3 = — £ Pi,P2,P3, e Ps) , B29) 

mt mt mf 

where the short notation £(pi,P2,P3, e(p3)) = £^^^"PiiiP2uPzp^a{pz)- The values of the coeffi- 
cients ^4] and A^- are 

A\ = (o?2(riV2 + ri^(2(w + x) + r2 + 4r2^) - r2(32; - 2v + r2{u -w + 2)-x) 

+ri(-i; + r2(-2 + 2u + ?,v + Ax + r2(-2 + 3r2) + ra^))) + dir4(2x(ri + 2r2) + 
+ r2 - rir2 + W)rA)) (B30) 



^1 ^ _2Limt^^ _ 2^^^^-^ ^ ^^-j ^ ^^^^ + n + ri^ + 3rir2 - r2(3 + 2r2))r4 + 



0^2(1 + r2)(M + (l + ri)r4)) (B31) 

7 

Al = (c^2(l + r,) + d,r,) (B32) 

9 

= -^^^(dir4(-2x(i; + ri + r2(2 + ra)) + {v{rx - 3r2) + r2(ri + ri^ + 3rir2 - 

2L2r4 2 ^ 

r2(3 + 2r2)))r4) + d2(4xVi + 2x(ri^ + r2(-l + w + r2) + ri2(-l + 3r2) + 

ri(-l + 2k + -y + 2r2^)) + r4((l + ri)r-2{2u + 2ri^ + (-1 + r2)r2 + ri(-l + 3r2)) + 

v{2u + {-l + r2)u)))) (B33) 

9 

^3, ^ ((ci2(2M(ii + x) + ri^ - (ii + 2a;)r2 - r2^ + ri^(-l + 3r2) + r2^{u - rg^) + 

2iv2vT4 ^ 
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ri^(2(u + x) + 3(-l + r2)r2 - r^^) + ri(-n - 2x + r2{3u + 2x + (-3 + r2)r2 - 
2r32))) + c/i(-4x^ + u{ri - 3r2)r^ + (-1 + n - 4r2)r4^ - 2x(m + 4:r2n)))) (B34) 



v4f = ^'^'' [{d2{2{u + x)+ 2n^ + 3nr2 + r2^ + r^) + di{-2x + {r, - 3r2)r4))) (B35) 



2L2^/rl 

And the values of the coefficients A? and Af are 



2 ^1 4L2mt2 / , 

= -— L2 3 — (2xct2'^2 - «2(-'i^ + (^^1 - 3r2)r2)r4 - dir2r4' 

Li r4 2 ^ 

= —^L2 ^-^(di{-2x + nn-3r2u) + d2{u + r4^^)) (B37) 



(B36) 



<9 4/72^2^*2 

= --^L2 (B38) 



v4f = -^L2 - ^^^^(rf2(riV2 + r2(2M - 2v + x + z) + n{z ~ 2v - 3x)) + 



di{ri - 3r2)r2r4 - 2xdi{r2 + r4)j (B39) 
= ^L2 - ^^(d,n' + d2r,{x + y)) (B40) 

= 4i^L2 + (B41) 

Here for convenience an overall factor Cg has been contracted out from these coefficients. 
Then the square of the amplitude |Mjp can be conveniently obtained with the help of 
Eqs. flB26]lB27llB28|) . 



2. Amplitude for spin-triplet S-wave state: P-Siji 

There are eleven basic Lorentz structures Bj for the case of (6c)pS'i], which are 

Bi = e(s^) ■e(p3), B2 = —e{pi,p2,e{s^),e{p3)), B3 = —e{pi,p3,e{s^),€{p3)), 

Bi = —e{p2,P3,e{s,),e{p3)), B5 = ^ ' ^6 = 2 ' 

o P2 ■ e{p3)p2 ■ e(s^) P2 ■ e{p3)p3 ■ e{s,) ipi ■ e{p3) 

Bj = 2 '^8= 2 '^9 = 4 e{pi,p2,P3,e{Sz>>, 

mf mf irq 

Bio = — £^(Pi,P2,P3,e(P3))P3 ■ e(-s^), Bu = —e{pi,p2,P3,e{p3))p2-e{s,), (B42) 
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where the polarization vector e{sz) related to the spin angular momentum of the spin-triplet 
state. The values of the coefficients Aj and A^' are 

Al = ^-^^(-{din{x-ri{z + r2)+r2{y + z + r2))) + d2{2u{v + x)+x{n-l) 

+2x{ri — l)r4 + v{2ri — l)r4 + 2r2^{u + (ri — l)r4) + r2{u — x + {u + x)ri + 

{2ri + w-2)r4 + (-l + n)n^))) (B43) 

7 

Al = ^^^^^(0^2(1 + 2m + n + r2 + nrs + 2r3^) + (ii(l + ri)r4) (B44) 



. '^4 

^3 = ( - rfir2r4 + (i2(4(v + + r2(l - ri + 2r2 + r^))) (B45) 

^4 = " (2ud2 + rfi(-ri + r2)r4 + rf2r4(-l - + 2r4)) (B46) 

7 

^1 ^ Z^^2^((i2(l + ri)(l+r2) + rfi(l + ri + 2r2)r4) (B47) 

7 

A\ = ( - dir2U + c?2(2(m + + r2(l - n + 2r2 + r4))) (B48) 

= -4diLimt^V^(l+r2) (B49) 

Ag = 2Limt^ ^/rl{d2 + diri - dir2 + d2r2) (B50) 

= zi^H^ (B51) 

Al, = ' i-* (B52) 



Af = //^^^^/_/ — (u + x)r2^ -\-{v - x)r^ + r2{x + + a;)ri + 

2Lo./r^^ 



2L2^/rl 

(1 + ri)r42)) + d2{2x{x + (u + :r)ri) + ,7;(2(t- - 1) - 2ri + r3^)r4 -{v-x){l + n)ri^ 
+vrA - r2^ri{u + r4 + riri) + r2{u + r4 + rir4)(2(M + x) + r4(-l + 2r4))))j (B53) 



^2 = ■7^f^^(-{dirA{2x + r2 + r2^ -ri{l+r2)+ri^))+d2{-2xri + 2{u + x + ur2) + 
2L2yJr 4 ^ 

r4(r3' + r4-rir4))) (B54) 

9 

A( = ;^^!^^((rfi(2a;-rir2 + r2^)r4 + rf2(2x(l+ri) -2Mr2 + 2t;r4 + 

r2(l+r2-2r4)r4))) (B55) 

Al = ^^tV^4 /^^^^2 ^ ^^^2(m + x) + r4(l + ri + 2r2 + r^))) (B56) 
2ivo ^ ' 
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9 

= — — ^—(dir4{2v - 2x + ri + r2 + rir2 + - r^^) + d2{2{u + x) - 2xri + 

r4(r3' + r4-rir4))) (B57) 
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Ag' = _,/^^^ ((c?i(2x - TxT2 + r2^)r4 + d2{r2r^ + 2(x + xri + vrA + 
2L2.y/ri^ 

r2(-2(M + x) +r4-2r4')))) (B58) 

= - ^^^"f^ (z. + 2x + r4^) (B59) 

^3' _ ^^^V^4 ^^^^^2 ^ ^^^2(m + x) + r4(l + ri + r4))) (B60) 

. *i (BCD 



= (B62) 
and A}y^ — ^l^Q = 0. And the values of the coefficients ^4^ and A^- are 

= -^^2 + ^^^^fdi(x - rir2 + r2^)r4 + d2{{v + 2x)r4 + 2r2V4 - 





ri(x + 2r2r4) + r2{y — 2u))^ 


(B63) 


42 - 

^2 — 


Al ^L2mfi , . N 
-— L2 ^^(rfir4 + ot2r4) 


(B64) 


42 - 
^3 — 


A\ Ad2L2mtir2 

2^2 ~r 


(B65) 


42 - 
^4 — 


A^ 7 

— -f^L2 + Ad2L2mt^y/rl 


(B66) 


- 

^5 — 


JL2 ^ (cZir4 + ^2^4) 


(B67) 




Al = 


Ad2L2mt^r2 

-t-/2 ~r 

Li ^/u 


(B68) 


A? = 


Al 7 
-—L2 - 8dlL2mt^^/U: 


(B69) 




j^i ^ 
-y^-^2 + 4(i2-^^2^t^"\/rl 


(B70) 



Af = 4^L2 - - 2x^d2 - xr4{-{di + 2^2) (ri - r2) + d2r4) + 

r4 (-1*^2^2 + vd2r4 + r2r4 i-din - d2T4))) (B71) 
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A 3' 

Ai = -^L2 - ^^(cii(ri - r2)r4 + ^2(2(1. + x) + r4')) (B72) 
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— — —{2x + r2ri) 


At 


Li 


, 9 3 

Ka2mt 2r4 2 


At 


43' 
Li 


9 

KITLt^ / , 2 7 / \\ 

- — —[diT^^ +d2{u + x + y)) 


At 


43' 


Kd2mfi 

— — —[2x + r2rA 


Ai 


Al 


Kd2mAJ 



(B73) 
(B74) 
(B75) 
(B76) 
(B77) 



and for j = 9, 10, 11 and m = 7, 9, 10, 11, we have 



Al A3' 

1? = -j^L2 and Ai = ^ 



Here for convenience an overall factor Cs has been contracted out from these coefficients. 
Then the square of the amplitude |Mjp can be conveniently obtained with the help of 
Eqs.f lB26]lB27llB28l) and Eq.(|33l). 



3. Amplitude for spin-singlet P-wave state: ^Pi]i 

The basic structures are similar to the case of ^5*1, only one needs to replace the po- 
larization vector e{sz) related to the spin angular momentum of the spin-triplet S-state to 
the present e{lz) that is related to the radial angular momentum of the spin-singlet P-wave 
state. Secondly, we need to change coefficients there to the present case. The values of the 
coefficients A] and A^- are 

L t 

A\ = — — ^—(-{di{ri-r2)ri'^{x + ri{r2 + z)+r2{r2 + y + z)))-2dr,rir2{2uri + 
rir2r4 2 

(r^ + rl- l)r4)(x + 2{v + x)r4 + 2r2^r4 + r2(2r4 + y)) + d2ri{ri'{x + yr2) 
+ri{x — 2u{v -|- x) -|- 2r2^ — r2^(— 2 + ?)U + v + r^^) + r2{v — u + ?)X + r^)) 
—ri{r-2{x — {w — 2)r4 -|- 2r4^) -|- ri{v + 2vr4 + 2x{l + Vi)) 

+r2\u + X + n{2 + 3n))))) (B78) 

Al = {-diin - l)(ri - r2)r4^ + 2rf5rir2(l + r2)i2un + {n^ + r^^ - I^a) - 

rir2r4 2 

c?2r4(-ri^(l + r2) + rir2(l + r2) + r4 + (r2 + 2{u + rs'^))^)) (B79) 
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■{d2ri{ri^r2 + {r2 + v + x + 3^)r4 - rir2(r2 + r^) + r2(o?i(ri - r2)r4^ + 



rir2r4 2 

2d5rir2(2wri + {rl + r| - l)r4))) (B80) 
im,2^r4 ^^^^^2 _ ^ ^^^^^^ ^ _ ^^^^ ^ ^ + r4'))) (B81) 

_Lirnt^^^^^^ + r2)(r2 + ri(l - n + r2))r4 + + n - 2r2)(ri - r2)r4^ + 

2rir2((i5(l + r2){2uri + {rl + rl - l)r4) + 2{d4{ri^r2 + r^(r2 - 2{v + x - 2z)) 
-r2(r2(u -w + 2)-2v-x + 3z)+ ri(-v + r2(-2 + 2u + 3v + Ax + 
r2(-2 + 3r2) + rl))) + dsr4{2x{ri + 2r2) + {v + r2 - rir2 + W)r4)))) (B82) 

3 {d2rA{ri'r2 + (r2 - - x + 32;)r4 - rir2(r2 + r^)) 
rir2r4 2 

+r2((ii(ri — r2)r4^ + 2d^ri{2v{2ri — l)r4 + 4(ri — l)r\ri + 4a;(ri + (ri — l)r4) 
+r2(2(M + 2x)ri + (-3 + 2^; - (-4 + ri)ri + r3^)r4 + 2(-l + ri)r4^)))) (B83) 
2Limt _ ^^^^"1^ +r2)r4 + 2rir2(d4(l + r2)(w + r4 + rir^) - d3{2x - 2uri + 



rir2^r4 

2r2V4 + r4^ - 2r4^ + n{u - 2(1 + r2)r4 + r^^) + r2(2r4 -u + 2y)))) (B84) 
(ci2(ri - r2)(l + r2 - r4)r4 + 4d5rir2(l + r2)(u + r4 + rir^)) (B85) 

(B86) 



rir2^/r4 

5 

-2d2Limt'^ Jrl 



rir2 
2Limt2 

rir2v/^ 
4Limt 



—unit 



24rir2(l + r2) + ci2r4) (B87) 

(o?4(l + r2) + ci3r4) (B88) 

7 

— 3- ((^1(^1 ~ r2)r4^(2x(K + x) — xri + {u + x)r2' + (x — v)ri^ + r2(x(ri — 1) 
4L2rir2r42 ^ 

+uri + (ri - l)r4^)) - 2d^rir2{2uri + (-1 + ri^ + r3^)r4) (-2a;^ + r4(-2Mr2 + 

+r2(l + r2 - 2r4)r4) - x{2u + r4(-2 + 2r2 + r4))) + (i2r4(2a;^(ri(-ri + r2) + r4) + 

x(-ri^(2w + n^) + r4(r2(2(-l + w + -u) + 2r2 + r^^) + r4^) + ri(2iir2 - (-2 + 2v + 

2r2 + r3^)r4 + (-2 + 3r2)r4^)) + r4(-vr4(-(ri - r2)(ri - r^) + r4) + 

r2{u{2u + r2 - r2^) + ri^(l + r2)r4 + (-1 + 2u + r2)r^ + uri{-l + r2 + 2r4) - 

rMr2 + r22 - 2(-l + r4)r4))))) (B89) 
7 

^ — j{2d^rir2{2uri + (-1 + ri^ + r3^)r4)(2a; + r^) + di{ri - r2)rA{2x + ri + 

4L2rir2r4 2 ^ 
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r2 + rir2 + r-i' + r^) + d2rA{2x{ri' — rir2 + r^) + r4(2ii(l + + r2r3^ + 

riV4 + - ri(r3' + r2r4)))) (B90) 

7 

T\di{ri — r2){2x + r2(ri + r2))r4^ + Axd5rir2{2uri + {rl + r| - l)r4) 

4L2rir2r4 2 

+d2r4(2w(ri - r2)r4 - 2x{ri{r2 - ri) + r^) + r2r4(ri(l + r2) - r2(l + r2) + 



2(« + r4')))j (B91) 

-/tmt^v^ /^^ _ ^^2 ^ 2d5rir2(2wri + (-1 + + r3')r4) + 
4L2rir2 ^ 

d2{-2{u + x)r2 + riV4 - (1 + 3r2)r4^ + ri{2{u + x) - r2r4 + n^))) (B92) 

7 

jidiiri - r2)(2a; - 2v + ri + ri^ - r2 + 3rir2)r4^ + d2ri{(2u + ri^ + ri^ + 

4L2rir2r4 2 

r2(r2 + ra^) - ri((-2 + r2)r2 + r3^))r4 + 2x{ri^ - rir2 + r4)) - 2rir2{d5{2uri + 
(-1 + ri^ + r3^)r4){u + x + y) + 2{d3r4{-2x{v + n + r2(2 + r2)) + {v{ri - 3r2) + 
r2(ri + ri^ + 3rir2 - r2(3 + 2r2)))r4) + d4{Ax'^ri + 2x{r/' + r2(-l + + r2) + 
ri^(3r2 - 1) + ri(2M - 1 + t; + 2r2^)) + + ri)r2{2u + 2n^ + (r2 - l)r2 + 
ri(-l + 3r2)) + v{2u + (-1 + r2)r4)))))) (B93) 



^{di{r\ - r2){2x + r2(ri + r2))r4^ + Adr,rir2{^x'^ri - (1 + ri)r2r4(-2M + 
4L2rir2r4 2 ^ 

(1 + r2 - 2r4)r4) + 't;r4(2w - (1 + ri)r4 + r4^) + 2x(2?/ri + (t; + (1 + ri)(r2 - l))r4 

+rir4^)) + d2ri{2v{ri - r2)r4 + 2x{ri^ - rir2 + (-1 + r2)r4) + r2r4(ri(l + r2) - 



r2(l+r2) + 2(ii + r4')))) (B94) 
:(di{ri - r2)ri{x + y) + 2dirir2{2u{u + x) + {u + 2x){ri - l)r4 + 



7 



2L2rir2y/rl 

uu^ - rz^ri + (ri - l)r4^) + 2d3rir2(-4x^ + 2x{-u + 2rir4 - 2r4(r2 + r^)) + 
r4(r4(r2 + r2^ - 4r2r4 - r4(2 + n) + ri(l + r2 + 2r4)) - 'w(2r2 + r4 - 2ri)))) (B95) 

7 

: (o?i(ri - r2)r4^ + 2(i5rir2(4ii('u + x) + {Ax{ri - 1) + 'u(4ri - 2))r4 



+{2u - 1 + ri^ - r3^)r42 + 2(-l + ri)r4^) + d2ri{-2{u + x)r2 + riV - 

(1 + r2)r4^ + ri{2{u + x) - r2ri + r^^))) (B96) 



^-d2r4' + 24ri(2iiri + (-l+ri2 + r3>4)) (B97) 



2L2rir4 2 
nd^rrit^ 
i^2s/rl 



(2(m + x) +r4(l + ri + r4)) (B98) 
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A'(^ = — ^—=(di{ri-r2)r^ + 2rir2{-2xd3 + ds{2ri-2r2-r4)r4 + d4{2{u + x) 

2L2rir2yrn^ 

+r4{l + n + u)))) (B99) 
And the values of the coefficients and Aj' are 

Al = -^L2 + '^^^^^^^^^{-di{ri-r2)r4^{x + r2r4) -2d5rir2{x + 2r2r4){2uri + 
Li rir2r4 2 

(-1 +ri^ + r^^)r^) + d2r^{xri' + r2r4(y — 2u — v^-2z)— ri{x{r2 + 2r4) 

+r4(i; + 2r2(r2 + r4))))) (BlOO) 

2^ 777- ^ 

Al = --^^2 + ^ *3 (4^d5riV2 + r4(ri(d2(ri - r2) + 24r2(-l + + rg^)) + 

ivi rir2r4 2 

(din- (rfi + d2)r2)r4)) (BlOl) 

A^ = -^L2 + ^^^H^^H!!!^ (B102) 
Li r-i 

Al = -^L, + ^L.rnfi^. _ 

Li rir2 

= --j^L2 3-(8mc?5'^i^''2 + 8d4rir2{-2xr2 + {-v + (ri - 3r2)r2)r4) + 

-t^i rir2r4 2 

r4(44rir2(-l + ri^ + rg^) + 2(-(diri) + dira + 4.d^nr2^)u + rf2(-2ri2 + 
3rir2 + r2^ + r4 + (-l + r2)r4))) (B104) 
Al = -^L2 + ^^^^(ci2r4' + 4ci5ri(2rix + r2r4) - r4(r2r4 -v + 2z))) (B105) 

A? = -4^L2 + ^^'"^!l ((ii(ri - r2)r4 + 2rir2(-2xci3 + 4(2ri - 2r2 - r4)r4 + 
Li rxr2s/rl 

d^iu + r^))) (B106) 
= -^L2 - '^Jf^{d2{r^-r2)r,-rU^r^T2{u-rT^)) (B107) 



= 4^L2 - (24rir2(2a: + r2r4)(2Mri + (-1 + n^ + r3^)r4) + 

Li 2Lirir2V^ 

di(ri - r2)r/i{xr\ + r2(y - -u)) + (i2r4(-2x^ - ri^r2r4 + vr^ — xr^' + 

2a;ri(r2 + r^) - r2^{y + x) + rir2{u + r4(r2 + 2r4)))) (BllO) 

< = ^i^2 - "^^P^idM - r2)u + rf2(y + U + X)) (Bill) 
ivi zLirir2 
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4 = - -J^Ill^(-2. + (r, - r,)r,) (B112) 

Li ZLirir2 

Li 2Lirir2 

^5 = —L2-\ 3-(2d4r2(-4a;ri(ri - r2)r4 + rir4(-2t;r4 + r2(4M + 3ri + 

-t^i 2Lirir2r4 2 

^2(2r2 - 1) + ri(5r2 - 1 - r4) + r^))) + r4(-x(8(i3rir2(ri + 2r2) + 2d2r4) - 
n{d,{n - f - 44ri(ri - 3r2)r2' + ^2(2^ + n^)))) (B114) 
A-t = - „ r _ (c?2(-2x + (ri - r2)r2)r4 - ^d^rxr2{vr^^ + r2V4 + 



Li 2Lirir2y^ 

ri(2x + r2r4) - 2?/r2)) (B115) 

4' = + ^^^;^v^ (4,,^ + d,{x + 2/)) (B116) 

At = - ^^p^id^n' + Ad,rM^ + y)) (BUT) 

< = ^L2 + (B118) 
Li Li 

At, = 4^L2 - '"^^7^^ (B119) 

= -4i^L2 and At, = ^L^. (B120) 

Here for convenience an overall factor Cg has been contracted out from these coefficients. 
Then the square of the amplitude |Mjp can be conveniently obtained with the help of 



and 



Eqs.( lB26]lB27llB28|l and Eq.(l33]). 



4. Amplitude for spin-triplet P-wave state: [^Pj]i 

There are totally thirty independent basic Lorentz structures Bj for the case of 
bc)[{^Pj)i], which are 

111 

Bi = — P2aep{p^)eip, B2 = — P2 ■ eip3)ei^, B3 = — pi ■ t{p3)ei^, 
rrit rrit rrit 

^7 = ^^£^(Pi,P2,a,e(P3))P3/3, Bs = ^^e{pi,p2,a,e{p3))p2f3, 
rrit rrit 

B9 = ^-^e{pi,p2,a,(3)pi- e{p3), Biq = ^^e{pi,p3,a,e{p3))p2i3, 



26 



Bi3 = ^-^£{Pi,P3,a,(3)p2- e{ps), Bu ^ ^-^e{p2,P3,a,e{p3))p2f3, 
nil 

-B15 = ^£(P2,P3,Q;,e(P3))P3/3, Bie^^-^e(p2,P3,a,P)pi- e(p3), 
mi rrii 

IS'^ ZS'^ ZS"^ 

Bi7 = — ^£(pi,Q;,e(p3),/3), Bis ^ —^e{p2,a,e{p3), P), -B19 = — ^£(p3, a, e(p3), /3) 

mf ml mi 

B20 = -^e{pi,p2,p3,a)pi ■ e(p3)p2/3, -B21 = -^e{pi,p2,p3,a)pi ■ e(p3)p3^, 
mi mi 

B22 ^ -fs{pi,p2,p3,-e{ps))psaP2,, B23 = ^e{p„p2,p3,-e{ps))psaP3,, 
mi mi 

= !ff.(p,.p2,p3..0,3)W... fta = ^P.-.feW3„ 
mi mi 

-B29 = -^P2 • e(P3)P3aP2/3, -B30 = -^P2 ' e(p3)P3aP3/3- (B121) 

2 1 

Noting the fact that e^p is the symmetric tensor and e^^ is the anti-symmetric tensor, 
and the fact that = e'^^ = 0, one may observe that the terms involving the following 
coefficients do not have any contributions to the square of the amplitude, so one can safely 
set the following coefficients to zero: 

(B122) 
(B123) 
(B124) 

The coefficients Aj, Aj, Af and A'j' that are the same for all the three [^Pj]i states (with 
J =1,2, 3): 

5 



A.i'Po) = 


for 


. = (1 


-4),j = 


(9,12,13,16,17,18,19) 


AjCPi) = 


for 


i = (l 


-4),J = 


(2,3,6,24,26,27,28,30) 


4('P2) = 


for 


i = (l 


-4),j = 


(2,3,6,9,12,13,16,17,18,19) 



^20 


Sd^Limi^ 


A^ 

^21 


A^ 

^23 


Sd^Limt^ 
\frl 


A^' 

^20 


A^' 
^22 


= A%, = ^ 





r4 

7 7 

-2Kd4mi^r2 ^3/ 2Kd^mi'^r2 

L2y/rA ' -^^2\/r4 



(B126) 
(B127) 



and for j = 20, 21, 22, 23: = -^^La and Af = ^Ls. 
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The coefficients A], Aj, Af and Af tliat are tfie same for botfi P-Po]i and [^Po\2' 

A] = ^^(cisCl + r2)(2uri + (r? + - l)r4) - d^rJx - ri(z + r2) + 

\/r4 



T2{t2 + ?/ + 2;)) + di{2u{v + x) + x{ri - 1) + 2x(ri - l)r4 + v{2ri - l)r4 + 
2r2^{u + (ri - l)r4) + r2{u -x + {u + x)ri + {w-2 + 2ri)r4 + (ri - I)r4^)pl28) 

A\ = zi^5^i!!!il(_a; + 2ii(T; + x) + {u - v - ?>x)r2 + (-1 + 3ii + v)r2^ - 3r2^ + 

ri^(r2 - V - a; + 32;) + ri(-v - x + r2(-l + 3v + 3a; + r2(-2 + 3r2)))) (B129) 
= (B130) 

A] = Zl^5^l!!^(i ^2u + nT2 + 2r32 + T,) (B131) 

^1 = ^b^^dz^lj^ri)ri + di{l + 2u + nr2 + 2r^^ + ri)) (B132) 



= -^^^^^{d^{r2 + v + x + 2,z)-dzr2rA) (B133) 

V^4 
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= l^^i^(,, + ^ + ^ + 3.) (B134) 
^1 = 1^2^(2^d4 + (rf3(-ri + r2) + d4(-l + 2ri + r2))r4) (B135) 



^1 ^ (^5(1 + ri)(l + r2) + d4(r2 - ^ - x + 3z) - d^r2n) (B137) 



^29 = -4Limt2y^((i4 + (ign - (i3r2 + (i4?^2) (B138) 
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= ^^^(d4(2x^ + 2^iV2 + 2riV2 + {u-v + ?,x)t2^ + (-! + « + 'i;)r2^ + r2^ + 
ri^(x + r2 + 5r2^) + ri^(— v — x + r2(3r2 -\- Au — 2>v + Az — 1)) + 
xr2(ii + + w - 2) + ri((ii - 2^; + 2x)r2 + (-2 + bu + 2v + 3x)r2^ + 
3r2^ + r2^ + x(m + V + w + 2a; - 2))) + (i3(-2a;(M + + :r/ 2 + /■i'"'/-2 + 

— ti — 2a;)r2^ + r2^ — a;ri(l + r2) + ri^(v — x + r2 + 2r2^) + rir2{u + 2v + 
r2{2 + r2)))r4 + 4(2iiri + (-1 + + r3')r4)(x + y)) (B139) 

= _'lp!Ihl{2x^ + 2?iV2 + (-1 + 2u + 2v)xr2 + 2riV2 + (m - ^; + 3a;)r2^ + 
(-1 + + 't;)r2^ + r2^ + riV2(l + 5r2) + ri^(-'t; - x + r2(-l + Au + v + 
3x + r2(3 + 4r2))) + ri(2x2 + + 2^; + r2(2 + 3r2)) + 
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r2{u -2v + r2(-2 + 5x^ + 2^; + r2(3 + ra))))) (B140) 



7 



At = ^^^(2«ri + (-1 + + r^y,) (B141) 

= _^^5!^(-2(m + x) + + ri^(-l + 2r2) + ri{2x + (-2 + r2)r2 - rg^) - 

r2(2ii + r2 + ra')) (B142) 

7 



= ^^^(rf4(-2(M + x) + ri^ + ri2(-l + 2r2) + ri(2x + (-2 + r2)r2 - rg^) - 

r2{2u + r2 + ra^)) + 4(20; + + ri(-l + r2) + r2 + 2r2^)r4) (B143) 

7 

= ^!!^(ci3(2xri + 2xr2 - riV2 + r2^) + ^4(2^ - 2riV2 + ri{2{v + x) + r2 - 3r2^) 
+r2(-2M + 2v + r2-r2^))) (B144) 
Afi = ^^!^(_2x + 2riV2 + r2(2^i -2v-r2 + ra') + ri(-2(^; + x) + r2(-l + 3<^^)15) 

.3/ Kmt5y^(o?3r42 + (i4(2(M + x) + 2ri2 + 5rir2 + 3r2^ + r4)) /d-, 
A4 = 7 (B14bj 

Af, = ^d^^t'^V^^ ^2{u + x) + {l + 2r, + 3r2)n) (B147) 
^2 
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= - ^^^^^ {d^{-2x + 2riV2 + r2(2(w - v + x) + {-1 + r2)r2) + ri{-2{v + x) + 
L2y/r2 

r2(-l + 3r2))) + d5{-2{u + x) + + rx^{2r2 - 1) + rx{2x + (ra - 2)r2 - r^) 
-r2(r2 + ri)) + (i3(-2x + (ri - r2)r2)r4) (B148) 
^3, ^ _ ^^*^v^ (^3^^2 ^ ^^(2(ii + x) + 2ri2 + 3rir2 + r2' + r4)) (B149) 



and 



A\ = -^L2 - ^^^^'' {d3{x - nr2 + r2yi - d^{2ur, + {r-,^ + rg' - 1)^) + 

di{{v + 2a;)r4 + 2r2V4 - ri(x + 2r2r4) + r2(t/ - 2u))) (B150) 



^4 = -^i^2 + !^^^^((^; + 2x)r4 + 2r2V4 - ri{x + 2r2r4) + r2{y - 2u)) (B151) 



= -^i^2 (B152) 

A^ 6 

= _- ^L2 -8d5i^2mt2 0^ (B153) 

= -^^2 + ^^^{dsr, + d,r,) (B154) 



,2 ^10 r Sd4,L2mt2r2 
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^11 — 


8d5L2TOt2r2 

1^9 + 

r /FT 

^1 5 

Li 


(B156) 
(B157) 


^15 — 


5 

— ^L2 + %d^L2mt^y/rl 

1 


(B158) 


^29 — 


^25 r SLsmtt 
— — -L2 ^— (d4r2 + d5'^4) 

— P^-^2 - 8d4L2mt^y/rl 
Li 


(B159) 
(B160) 


Af = 


-y-L2 - —^{2x^d4 + xn{-{{d3 + 2(i4)(ri - r2)) + (i4r4) + 
Li LiWr4 





r4(-('i;d4r4) + dzr2ri + d4r2{u + r^))) (B161) 

Aj = -7^1/2 — ^(a;(2ri - 2r2 - r4)r4 - 2a;^ + r4(t'r4 - tir2 - r2r4^)) (B162) 

At - ^L2 (B163) 

= ^L2 - ^^^(2(« + + r4^) (B164) 

< = 4^L2 + ^^^(4(ri - r2)r4 + c^4(2(« + x) + r4')) (B165) 

< = #i^2 - ^^^(2a; + r2r4) (B166) 
Af, = ^L2 + ^!l^(2x + r2r4) (B167) 
Al = 4liL, + ^J^k^ (B168) 

^25 = 4^^2 - ^^^(2x^4 + 2x(d4 + 4) + d4r2U + 4^4') (B170) 
Li Liy/r4 

Li Li 

3, 



The remaining coefficients for the case of [ -Poji^ 



5 



A\ = (^2(1 + r2)(ri^ - ri^ + r2(^ + ra) + ri{-u + r2(2 + ra))) + rfi(r/ + 

2xri(l + ra) + ri^(7r2 - 1) + rir2{6u + (r2 - 7)r2) + ri^{u + r2(9r2 - 5)) - 
r2(2x(l + r2) + r2(r2(3 + 2r2) - u))) + 2d5rir2(l + r2){-2uri - 
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(ri'-l + r3>4)) (B172) 

— ^^^—^{d2{-ri'r2 + ri^(r2 + v + x + z) + rir2{2x + r2(2(M + v + 2x) + 
r2(3r2 - 1))) + ri^{v + r2(3r2 — u — v — w — 3x + 5z + 2) — r2'^{v + 2x + 
3r2^ - r2(-2 + v + rs^))) + c?i(riV2 + r2^(r2 - 3v + Az) + ri^(v + 2x + 
r2 + 6r2^) + rir2(6r2 — Tt" — 3a; + 132;))r4 + 2(i5rir2(2(t' + x) + 

r2{2 + ri + 3r2))(2iiri + (-1 + + r3')r4)) (B173) 

L mi 

^^(-4(ri - r2)(l + r2) + 0^1(^2 + Qnr2 + r2')) (B174) 
rir2V^ 

(B175) 

^^^^^^*-i(r2-.-x + 3.) (B176) 

(o?4(l + ri) (1 + rs) + (i3(l + ri + 2r2)r4) (B177) 

/r4 

8ci3LimttV^(l + r2) (B178) 

44ivimtiV^(l + r2) (B179) 
7 

— ^^(2(i5rir2(2tiri + (ri^ - 1 + r3^)r4)(a; + + (ii(4a;^(r2 - ri) 

4-L2rir2v^ 

+ ri^(l + 6r2) + r2(3ii + r2 + 4r2^) +ri{u + r2(6 + 13r2)))r4^ + 

2a;(M(r2 — ri) + 4r2r4^)) + d2{r2{—Au{u + x) — {u + 2x){2ri — l)r4 — 2Mr4^ + 

2r3W + (1 - 2ri)rA^) + ri{-r^^ri^ + ri{l + rA){x + y) + u{u + x + y)))) (B180) 
7 

3 (2x(ii(ri' - ri(v + (-5 + r2)r2) + r2(v + r2(2 + r2)))r4 + 
4L2rir2r4 2 

di{v{ri + 3r2) + r2(ri - ri' - 5rir2 + r2(3 + 2r2)))r4^ + 2d;,rir2{2uri + 

(-1 + ri^ + r3^)r4)(2x + vn + r2V4 + r2(r4 - 2y)) + d2(4a;Vi(ri - r2) + 

2x{ri^ + ri^ - r2^(v - 1 + r2) + ri^(-l + 2m + v - 3(-l + r2)r2) + 

rir2{-2u + r2 - 2r2^)) + v{ri^{l + r2) + 2ri{u + r2) + 

r2(-2ii + r2 - r2^))r4 + r2r4(-2ii(r2 + + r4)) + r4(ri2 - 2^^ + 

2rir2 - 5riV2 - r2^ - rir2' + r^)))) (B181) 

-(/cmt2 ^^^^^ _ ^ d2{2ri^ + ri^{r2 - 1) - 2ri{-u -x + r2 + r2^) 
AL2nr2^A 

-r2(r2 + 2u-v + x + z))+ di{ri + 3r2)r4^ + 2d5rir2{2uri + (n^ - 1 + rs^li^l'^) 
2^Wv^ (B183) 
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_ Kd^mt^ (2(a; + {v + x)ri + ^2) + r2^ri + r2(-2(M + x) + r4 - 2r4^)) (B184) 

7 

-^^^(-(4(-2w + 2a; - r2 + ri(-l + r4))r4) + di{2{u + x) - 2xn + 

r4(r3' + n- rm))) (B185) 

'"^^7*^^ (0; + y) (B186) 

/^^rn^^yFI ^^^^ + x) + (1 + 2ri + r2)r4) (B187) 



^1 2L m ^ 

--^L2 ^^(-2xdi(-2r2 + r4) + 2d5rir2{2uri + (-1 + + r^y^) + 

J-'i rir2^r^ 

dxr^{2r2r^ + r^) + d2{-2ur2 + ^4 - r2^r4 + rir4(-3r2 + r^))) (B188) 

_^L2 + ^^^^(2xci2(ri - r4)(-2r2 + r4) + ci2r4(^;(2r2 - r^) + 
L\ rir2r4 2 

(ri - r2)r2(-2ri - 4r2 + r4)) + r2(-(dir4(2(ri - r2)r2 + (n + 3r2)r4)) + 



^d^rxr2{-2urx - (-1 + + r3^)r4))) 


(B189) 


2d2L2mt^ 

-—L2 — (-2r2 + ta) 

Li rxr2\/rl 


(B190) 


^24 T 


(B191) 


r , 8d5L2mt^r2 

1^2 + 


(B192) 


8L2mti 


(B193) 


-^L2 + \U^L2mt^^^ 


(B194) 


— -^L2 + 8d5L2mt^y^ 


(B195) 


^L2 + '^!^''^{dir42{n - r2)r2 + {n + ^r2)u) + d2n{x + y)) 


(B196) 



At 



1^L2 , {2d5nr2{2x + r2r4){2uri + (-1 + + rs^y^) + 

-^1 2Lirir2r4 2 

(iir4(r2r4(2r2r4 + r4^) + 2a;(2r2r4 + ri(2r2 + r4))) + d2r4{2xri'^ + 2xrir2 + 



nu{v + 3r2^ - r2r4) + r2(r2V4 - 4yr2 + r4{2u + v + 2x + 2r4^)))) (B197) 

^L2 - ^^^H^ (B198) 
Li 2Lirir2 
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^24 T 




^26 ~ 


43' 


r ^ (2x + r2r4) 


^27 ~ 




7 

- — ^(4^4 + ct4(ii + a:; + y)) 


A^' - 

^30 — 


43' 

^30 7- 





(B199) 
(B200) 
(B201) 
(B202) 



The remaining coefficients for the case of [^Pi\i state: 
A\ = ^—^{d2{l + r2){u + r2 + ri{-l + r4))r4-dir4^{u+{-l + n)u) + 



2rir2{—{d4{—x + 2u{v + x) + {u — v — 3x)r2 + ri^r2 — 3r2^ + ri^(2(w + x) + 
r2 + 4r2^) + r2^{2u + w - 2) + ri{-v -x + r2{u - 2r2 - v + w + 3z - 2) + 
rs'^W + dsix + riVa + ri{-v - x + (r2 - l)r2) + r2(r2 + u-v + 2z))r4))(B203) 

Al = '^^^^*_ {d2{-x + ri(r2 + v + x + z) - r2{r2 -v + 3^))r4 + x{2d^rir2{2u + 
rir2v^ 

(ri - 1)(1 + 2ri + 3r2)) + dW) + r2(dir4^ + 2d^n{ur2{l + 2ri + 3r2) + 

v{2u + ri(2ri - 1) - r2 + 3rir2 + r22) + r2((ri - 1)(2 + n + Sra) + r3^)r4)092O4) 
^ 2Lim,2^ _ ^^^^ ^ ^ ^ 

A] = ^^^i^^^(4Mci5ri(ri - r2)r2 + r4(2ci5ri(l + ri)(-3 + 2ri - r2)r2 + 

rir2r4 2 

d2(2ri^ - ri(-3 + r2) + r2(3 + r2)) + di(ri - r2)r4)) (B206) 
^8 = ^^^^'' {rir2{d^{l + ri)r4 + 0^4(1 + 2u + rir2 + 2r^ + r^) - + r2)r4)(B207) 

^l'"2V'"4 

^1 ^ Z^^2!!^(o?2(ri(2 + n) + r2)r4 + 2d^r^r2{2ur^ + (-1 + ri' + r3')r4)) (B208) 
r\r2r^ 2 

^lo = ^^^i^(-(i2(ri(-3 + r2) -3r2(l + r2))r4 + o?ir4^ + 2rir2(4W5ri + 
rir2r4 2 

2d5(-l + + r3^)r4 + r4(d4(r2 + t; + x + 3z) - d3r2r4))) (B209) 
A\, = ^^^^ir2 + v + x + 3z) (B210) 

"^^^ - ^^^^^^ 

^}3 = ^^(8W5riV2 + r4(d2(-(ri(-3 + r2)) + 3r2(l + r2)) + 
rir2r4 2 
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4d5rir2(-l + ri' + rs') + di{n - r2)r^)) (B212) 
A^^ = ^^^^^{2udi + (4(-ri + r2) + ci4(-l + 2ri + r2))r4) (B213) 
4L m i 

= ' ' (-2W5rir2 - (rf2(ri - r2) + d^rM-l + 2ri + r2))r4) (B214) 
rir2V^ 

A\q = i^tV^4 ^^^ _ ^^^^ ^ 2d2ri + (iir2 - 4r2) (B215) 
rir2 

= - 3 (c^2(-riV2 + ri^(4(t; + x) + r2) + r2(v + 4a; + 3r2^ - r2(-2 + 
rir2r4 2 

2m + 3v + ra^)) + ri(v + 4x — r2(M — v + 'u; + x + 3^ — 2)))r4 + 

di{n\2 - r2{v + r^) + ri(r2 - 2^; - x + 3^))r4^ + 2dr,rir2{2{v + 2x) + 

r2(2 + n + 3r2))(2Mri + (-1 + + r3^)r4)) (B216) 

A\s = ^''"*3 (-24rir2(-2M - (-1 + 2ri + r2)r4)(2iiri + (-1 + r^^ + r3')r4) + 
rir2r4 2 

dir4^(M(ri - r2) - 2x{l + r2) + (ri - 1 - 2r2)ri) + (i2'^4(M(4ri^ - ri(r2 - 3) 

-(^2 - 3)r2) + r4((ri - r2)(l + r2 + 2r^^) + ri(l + 2ri - r2)r4 + 2r4^))) (B217) 
L mi 

Aw = -^—^{d2ixri{-3 + r2) - 2rir2^ - ri^(r2 + f; + x + 3z) + r2(-3x + 

rir2r4 2 

r2(r2 + ?>v + z)))ri - dir^{x + r2r4) - 2dr,rir2{2x + r2r^{2uri + 

(-l + ri2 + r3>4)) (B218) 

^25 = ^^^^(c?2(l + r2)r4-ciir4' + 2rir2(4(l + ri)(l + r2) - 

(i4(r2 - - x + 3z) + o?3r2r4)) (B219) 

^29 = -4Limti V^(d4 + 4n - (^3^2 + rf4r2) (B220) 

7 

= — —{2dirir2{-2x{x + {u + x)ri)+x{2-2v + 2ri-r^^)ri + 

2L2T\T2Jrl 



{v — x){l + ri)r4^ - vr4^ + r2^r4{u + r4 + riVi) + r2{u + r4 + riVi) ■ 
(-2(u + x) + r4 - 2r4^)) + r4(2a;cii(-M + r4^) + dir^^iu + 2v + n + 

+ ?^i?^2 + ^^2^ + + 0^2 (2m^ + 2^0; + Mr4^ - r3^r4^ - r2(a; + + 
ri(l + r4)(a: + y)) - 2d^nr2{-{x{2{u + x) + ri)) - + x)r2^ + 

- x)r4^ + r2{x + {u + x)ri + (1 + ri)r4^)))) (B221) 

7 

^4' = ^ ^ (24rir2(2x(x + + x)ri) + x{-2 + 2v- 2ri + r3^)r4 - 

2L2rir2V^ 

(■u - a;)(l + ri)r4^ + vr^ — r2^ri{u + r4 + riTi) + r2(ii + r4 + rir4) • 
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{2{u + x)+ r4(-l + 2r4))) + r4(di(2x^ + xr^^ + (n - r2)r2r4^) + 

d2(-(x(2(M + x)+ ri)) + xrir4 + {2v - x)r4^ - r2^{2{u + x) + r^^) + 

r2{x + 2{u + a;)ri + (1 + 2ri)r4^)))) (B222) 



''' (^2(2(^ + x) + ri(2ri - 1) + r2 + 3rir2 + r2') + di(2x + r^)) (B223) 
zL2r ir 2 

7 

^-(dir4(2x + r4^)) + d2r4(2(ii + x) +r2 + ri(-l + r4) +r4^) 



2L2rir2^/rI 

-24rir2(2a;(ri - 1) + 2M(ri - 1 - r2) + (ri^ - 1)^ + (n - l)r4^)) (B224) 

7 

— ^^^^^^^^(2d3rir2(2x + + ri(r2 - 1) + r2 + 2r2')r4 - 2dir4(2x + r^^) 
ZL2r ir 2\/Ti^ 

—2d/^r\r2[^—2xr\ + 2{u + x + ur2) + ?"4(?"3^ + r4 — rir4))) (B225) 

7 

3 (24rir2(l + r2)(2uri + (ri^ - 1 + r^^)r^) - d2r^{-{r^^{\ + r2)) + 
2L2rir2r4 2 

(1 + r2 - r3^)r4 + ri(2x + r2 + r2^ + r^))) (B226) 

7 

— —iAxdirir2 + Axd4ri'r2 — Aud^riT'i' — 2xdir4 + Axdzr\r2rA + 

2L2rir2v^ 

Avd4rir2r4 + 2d4rir2r4 — 2d3ri^r2^r4 + 2d3rir2^r4 + 2d4rir2^r4 - 

Adirir2^rA^ + dir^" + 2d^rir2{2uri + (-1 + ri^ + r3^)r4) - 

d2ri{2{u + x) - r2(l + 2r2) + r^^ + ri(l + 2r2 + Vi))) (B227) 

7 

^((^2(^1 - r2)r4 + d^rii2x{l + ri) - 2Mr2 + 2wr4 + 

L2ri.Jrl 

r2(l + r2-2r4)r4)) (B228) 
7 



.Krf2(riV2 - 2ri(z + r2) - r2(2^; + r2(2 + r2)))r4 + 
2L2rir2r4 2 

2d^rir2^{2uri + (-1 + + r3^)r4)) (B229) 

~^'^'^'' {d2{-2{u + x) - ri(l + 2ri) + r2 - 5rir2 + r2')r4 + 
^L2r ir 2^r 4 

24rir2(2Mri + (-1 + ri^ + r3^)r4) + dir4(2a; + r4^))) (B230) 
_ /tm,2^ _^^^2 ^ ^^(2(ii + x) + 2ri2 + 5rir2 + 3r2' + r^)) - d^r^) (B231) 



Kd^mt-^ ^/^(2(x, + ^) + (1 + 2ri + 3r2)r4) (B232) 



2 



:(ci2(2(M + x) + 2r^ + ri(-l + r2) + r2 + 3r2^)r4 + 



4L2rir2^ 

2d5rir2(2^/ri + (-1 + + r3^)r4) + dir4(2x + r4^)) (B233) 
7 

— ^^^^^^^^ — j^{2xd\{v + ri + r2^)r4^ + + r2 - rir2 + 2r2^)r4'' + 
4L2rir2r4 2 
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2d5rir2i2uri + (-1 + ri^ + r3'^)n){Ax + n{v + + - 2r2n)) + 
d2ri{v{-2u + 3ri(l + rg) - r^iZ + r2))r^ - r2{2u + 2ri^ + ri^(r2 - 1) + 
r2 - ri(l + r2(4 + rs) + 2r32) + r2(-r2 + 2(m + r3^)))r4 + 

2x(2ri(ri - rs) + (-2 + + rsy^))) (B234) 



7 



* :(2o?5rir2(M + 2x + ri(2 + n) + 2r2 + 4rir2 + 3r2^)(2Mri + 



4-L2rir2A/rI 

(-1 + ri^ + r3^)r4) + d2ri{-2u^ - u{2ri^ + ri(-3 + rs) + 3r2(l + r2)) - 
2x(?i + r2 + ri(r4 - 1)) - ((n - 2)(ri - l)ri + 2r2 + 4rir2 + 3(1 + ri)r2^)r4 
+^3^4^) + dir^^Ax^ + (m + ri + ri^ - r2 + 3rir2)r4^ + 2x{u + 2r4^))) (B235) 

. (c?2(2x(^ + x) + xr2 + 4riV2 - {Au + 2^; + x)r2^ - 3r2^ - 3r2^ + 

4L2r ir 2^Jr 4 

ri^(— 2f + (—3 + r2)r2) — ri{x + r2(6r2 — Au — 2v — 9x + Qz)))ri + (ii(2x + 

ri' — r2^)r4(x + r2r4) + 2d^rir2{x + 2r2r4) {2uri + (ri^ - 1 + r3^)r4))) (B236) 

' -Aud4rir2^ + Avdirir2ri + 2c?4rir2^r4 - 2d3ri^r2^r4 + 



2L2rir2./rI 

2d^rir2^r4^ + 2(i4rir2V4 - 4(i4rir2^r4^ + dir^ + (i2r4(2(M + x) + r2 + 

ri(-l + r4) + r4^) + 2dr,{-l + ri)rir2(2(u + x) + ^{1 + n + r4)) + 

2a;(-rfir4 + 2rir2(o?4 + o?4ri - to + to)))) (B237) 
7 

- ^7*' (^i^2(to^ + ci4(2(K + x) + 2ri2 + 3rir2 + r2^ + r^) - dir^^)) (B238) 

-t>2^ 1^ 2 



Li rxr2Jrl 



rx{z + r2^) + r2(3z -u- 2v))) + d2{u + + r2^)r4 + d^ri") (B239) 

AL2mt^ , , , , . . 

-1-^2 + ——={d2{-x + (n - r2)r2)r4 + 

J^i rxr2^ri^ 

2(i5rir2(-(t; + 2a;)r4 - 2r2^r4 + r^ix + 2r2r4) - T2{-u + a; + r^))) (B240) 

Ail; ^ 4(i2-^>2^i^\/^ , X 

-12 H ^ B241) 

Li rir2 

- ^^'"X (-3(i2r4 + 2d^rxr2r^ (B242) 

_^L2 + ^^:H!!^(-(ciir4) + d^T^T2T^ + tor2r4) (B243) 
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.^L^ _ l^:^(2d,rir2^ + U^r,) (B245) 
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L2 (B247) 



Li rir2 



^1 



L2 - 8d4L2mt^^ (B249) 



41 5 

■-^L^ + Sd^L^mt^^ (B250) 



Li rir2 



5 



--j—L2 3-(«2r4(2ri r2 + + 4a; + 3r2 )r4 - rir2(2r2 + r4)) + 

r2(o?i(ri - r2)r4^ + Mzrir2{2uri + (-1 + ri^ + r^^y^))) (B252) 
^18^2 + ^^^"^*l (-2rf5rir2(2^ri + (-1 + + r3^)r4) - rfir4(2a; + n^) 



Li rir2y/rl 

-d2ri{3u - r2^ + 2r4^ + ri(r2 + r4))) (B253) 



Li rira 



5 



^L2 - ^^^(d2r4 + 2rir2(-(d4r2) + d,r,)) (B255) 



^29 - - — 5 



^1 



L2 - 8(^4-^2^(2^ (B256) 



1/2 + ^(2(i4rir2(-2a; - xr4(-2ri + 2r2 + r4) + r4(-iir2 + i;r4 

Li Lirir2^r4 

-^2^4^)) + r4,{dir4^ + d2(?'i - r2)(a; + + 2dsrir2{xri - r2{x + r4^)))) (B257) 

— L2 + -p^\-d2r/^[-xrx + r2(x + r4 )) + 

Li ivirir2VrI 

2rf5rir2(-2x^ + x(2ri - 2r2 - r4)r4 + r4{-ur2 + vu - r2r4^))) (B258) 

— L2 — (-ri + r2) (B259) 

Li Lirir2 

/1 3' - 

-^L2 - (d2(r2 - ri)r4 + 2d5rir2(2(« + x) + n^)) (B260) 

Li ivirir2V^ 



-^L2 + — ^((i3(ri - r2)r4 + d^iu + x + y)) (B261) 
Li 
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/1 3' - 

^9 = — 1/2 - 3:("2r4(ri -rir2-r2r4) + 

2d5rir2{2uri + (-1 + + r3^)r4)) (B262) 
44' _ ^10 r Kmti(d2(-ri +r2)r4 + 2d4rir2(2x + r2r4)) /Tao«Q^ 

Af, = 4^L2 + ^^^(20; + r2r4) (B264) 

Li Liri 

^3 = ^i^2-^^^^(n-r2) (B266) 

< = #L2 + (B267) 

Ai = 4&L2 - ^!^^^!^ (B268) 
Li Li 

Li 2Lirir2 

Afr = -F^^2 - 3 (2d5rir2(4x + r2r4)(2Mri + (-1 + + r^y^) + 

J^i 2Lirir2r4 2 

diri{2xri + r2r4^) + d2ri{Axri{ri - r2) + riri{'iv + r2^ - r2r4) - 

r2ri{Au + 2v + x-2y + z))) (B270) 
7 



^^8 = + ^?^^^^(rfi(ri - r2)r4^ + Ad.rM'^ur, + (-1 + + r32)r4) + 

ZLiT iT 2 



4' 
■18 


}'L2 + 

d2{—2xr2 


■19 — 


^19 T 


4' 

■25 — 


43' 

^25 T 


4' 

■29 ~ 


''L2 + 



2„ Q„ , ^ 2\ , ^ /Q„, ,0^ O I „ 2\ 



(B271) 

A% = ^L2 - !!^!^^ (^,, _ + 3r2r4^) (B272) 



7 

^25 = ^^2 - (d2(ri - r2)r4 + 2rir2(4(« + a: + - (i4(2x + r2r4))0B273) 

Lirir2V^ 

A- = :^L2 + (B274) 

Here for convenience an overall factor Cs has been contracted out from these coefficients. 
Then the square of the amplitude |Mjp can be conveniently obtained with the help of 
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Eqs. flB26]lB27llB28|) and Eqs.flMl [35]l36D . 
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